Fate of Zero-Temperature Ising Ferromagnets 
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We investigate the relaxation of homogeneous Ising ferromagnets on finite lattices with zero- 
temperature spin-flip dynamics. On the square lattice, a frozen two-stripe state is apparently reached 
approximately 1/4 of the time, while the ground state is reached otherwise. The asymptotic relax- 
ation is characterized by two distinct time scales, with the longer stemming from the influence of 
a long-lived diagonal stripe "defect". In greater than two dimensions, the probability to reach the 
ground state rapidly vanishes as the size increases and the system typically ends up wandering 
forever within an iso-energy set of stochastically "blinking" metastable states. 
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What happens when an Ising ferromagnet, with spins 
endowed with Glauber dynamics 0|, is suddenly cooled 
from high temperature to zero temperature (T = 0)? A 
first expectation is that the system should coarsen H] 
and eventually reach the ground state. However, even 
the simple Ising ferromagnet admits a large number of 
metastable states with respect to Glauber spin-flip dy- 
namics. Therefore at zero temperature the system could 
get stuck forever in one of these states. 

In this Letter, we present evidence that the behavior 
of such a kinetic Ising model is richer than either of these 
scenarios. While the ground state is always reached in 
one dimension, there appears to be a non-zero probabil- 
ity that the square lattice system freezes into a "stripe" 
phase, at least for equal initial concentrations of "f and 
I spins The relaxation is governed by two distinct 
time scales, the larger of which stems from a long-lived 
diagonal stripe "defect". On hypercubic lattices {d> 3), 
the probability to reach the ground state vanishes in the 
thermodynamic limit and the system ends up wandering 
forever on an iso-energy subset of connected metastable 
states. Again, the relaxation seems to be characterized 
by at least two time scales. It bears emphasizing that 
these long-time anomalies require that the limit T — > 
is taken before the thermodynamic limit L — s- oo' very 
different behavior occurs if L ^ oo before T ^ Q| . 

We can easily appreciate the peculiarities of zero- 
temperature dynamics for odd-coordinated lattices, such 
as the honeycomb lattice. Here any connected cluster in 
which each spin has at least 2 aligned neighbors is ener- 
getically stable in a sea of opposite spins. For any ini- 
tial state, a sufficiently large system will have many such 
stable defects. Because the number of such metastable 
states generally scales exponentially with the total num- 
ber of spins N , the system necessarily freezes into one of 
these states. However, on even-coordinated lattices the 
number of metastable states grows as a slower, stretched 
exponential function of N , and they affect the asymptotic 
relaxation in much more subtle way. 

We therefore study the homogeneous Ising model, with 
Hamiltonian TL = —J^i^ij-^ '^iO'j, where Ui = ±1 and the 
sum is over all nearest- neighbor pairs of sites {ij). We 



assume initially uncorrelated spins, with aj{t = 0) = ±1 
equiprobably, which evolve by zero-temperature Glauber 
dynamics m, corresponding to a quench from T = oo 
to T — 0. We focus on d-dimcnsional hypercubic lat- 
tices with linear size L and periodic boundary conditions. 
Most of our results continue to hold for free boundary 
conditions and on arbitrary even-coordinated lattices. 

Glauber dynamics at zero temperature involves pick- 
ing a spin at random and considering the energy change 
AE if the spin were flipped. If AE < (>0), the flip is 
accepted (rejected), while if AE = 0, the attempt is ac- 
cepted with probability 1/2. After each event, time is up- 
dated by l/L"^, so that each spin undergoes, on average, 
one update attempt in a single time unit. In practice, we 
update only the flippable spins (those with AE < 0) and 
update the time by 1/ (number of flippable spins). For 
each initial state, one realization of the dynamics is run 
until the flnal state. At T = 0, metastable states in this 
dynamics have an inflnite lifetime and these can prevent 
the equilibrium ground state from being reached. This is 
the basic reason why dynamics at T = is very different 
from that of small positive temperature. 

In one dimension, it is easy to determine the ultimate 
fate of the system [^. In T = Glauber kinetics, the 
expectation value of the i^^ spin, Si = {<Ji), obeys the 
diffusion equation and therefore the average magnetiza- 
tion (m) = i Ylfj is conserved |0] . Since there are no 
metastable states in one dimension, the only possible fi- 
nal states are all spins up or all spins down. For initial 
magnetization m(0), a flnal magnetization m(oo) = m(0) 
can be achieved only if a fraction ^(1 + m(0)) of all re- 
alizations of the dynamics ends with all spins up and a 
fraction i(l — m(0)) with all spins down. 

On the square lattice, there exists a huge number of 
metastable states which consist of alternating vertical (or 
horizontal) stripes whose widths are all > 2. These arise 
because in zero-temperature Glauber dynamics a straight 
boundary between up and down phases is stable; a rever- 
sal of any spin along the boundary increases its length 
and raises the energy. Note that a stripe of width one is 
not stable because it can be cut in two at no energy cost 
by flipping one of the spins in the stripe. 
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FIG. 1. Probability that an L x L square eventually reaches 
a stripe state as a function of 1/ logj L for L up to 1024. Each 
data point (with error bars) is based on 10^ initial spin con- 
figurations for L < 256, and then 13,500 and 2,700 configura- 
tions respectively for L = 512 and L — 1024. 

The mere existence of these metastable states implies 
that a finite sample may not reach the ground state. 
However, one might expect that the probability to reach 
the ground state approaches unity as the system size 
grows: lim^^oo ^gs = 1- Our numerical simulations on 
L X L squares appear to disagree with this expectation 
(Fig. |l|) . The probability Pgg grows extremely slowly with 
L and extrapolates to a value of approximately 3/4 as 
L — > cxD, so that the probability to reach a two-stripe 
state, Pstr = 1 — Pgs would be non-zero. We also find 
that states with more than two stripes almost never ap- 
pear. 



2.5 



1.0 
0.5 
0.0 



L=24 i 




(a) 




(b) 


■ L=50 fc 
L=100 / 


m 










i 

t 


1 

i 






°° oL=20 '\ 


i 






■ L=24 V 


't 

i 


M 
% 

\ 






L=30 ^ 


■ 




f 


— y=3(l-ml/4 \, 



-1.0 -0.5 0.0 
m 



0.5 



-0.5 




m 



0.5 



0.8 
0.6 
0.4 
0.2 




FIG. 2. Final magnetization distribution on the (a) square, 
and (b) cubic lattices. On the square lattice, this distribution 
narrows as L increases, but appears to reach a non-singular 
limit. On the cubic lattice, data collapse occurs even at small 
sizes. The number of realizations is 10^ for the square and 
> 10'' for the cubic lattices. 

When the two-stripe state is reached on the square 
lattice, both stripes have width typically of the order of 
L/2, as seen by a gradual narrowing of the continuous 
component C2d{Tn) of the final magnetization distribu- 
tion, F2d{m) = iPg, [S{m - 1) + S{m + 1)] + PstrC2d(m) 
(Fig. §(a)). However, as L increases the magnetization 
distribution appears to converge to a finite-width scal- 



ing limit. On the simple cubic lattice, there are many 
more metastable state topologies and also relatively more 
states with narrow stripes, so that there is a larger prob- 
ability that the final magnetization is close to ±1. The 
final magnetization distribution also exhibits good data 
collapse even at relatively small system sizes. Strikingly, 
the final magnetization distribution on the cubic lattice 
is well fit by F^dim) = |(1 - m^) (Fig. |(b)). 

Intriguing behavior is also exhibited by the survival 
probability S{t) that the system has not yet reached its 
final state by time t. As shown in Fig. |[ S{t) is controlled 
by two different time scales. On a semi-logarithmic plot, 
S{t) lies on one straight line with large negative slope for 
intermediate times and crosses over to another line with 
smaller negative slope at long times. In this intermediate 
time regime, the energy decays as as expected |^ 

The crossover in S{t) occurs when domains, which grow 
according to the classical t^^'^ law reach the system 
size. This gives the crossover time Tc (x L^. 

Quite surprisingly, the source of the long-time anomaly 
in S{t) arises from the approximately 4% of the configu- 
rations in which a diagonal stripe appears (Fig. ^). On 
the torus, this configuration consists of one stripe of t 
spins and another of | spins which, by symmetry, have 
width of order L. Each of these stripes winds once in 
both toroidally and poloidally on the torus; they cannot 
evolve into straight stripes by a continuous deformation 
of the boundaries. Consequently a diagonal stripe con- 
figuration should ultimately reach the ground state. 
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FIG. 3. Time dependence of the survival probability S{t) 
on Lx L squares. Main graph: S{t) versus t/Mio to highlight 
the long-time exponential tail, where Mk = (t'")'''* is the k*'^ 
reduced moment of the time to reach the final state. Scaling 
sets in after S{t) has decayed to approximately 0.04. Inset: 
S(t) versus t/Mi/iQ to highlight the scaling and the faster 
exponential decay in the intermediate-time regime. 

Diagonal stripes are also extremely long-lived. For 
L — 200, for example, the time for such a configuration 
to reach the final state is two orders of magnitude larger 
than the typical time. To understand this long lifetime, 
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we view a diagonal boundary as an evolving interface in 
a reference frame rotated by 45° |6) . 




FIG. 4. Diagonal stripe configuration on the square lattice. 
The lower portion shows part of one interface rotated by 45° . 
Zero-temperature Glauber dynamics allows for either particle 
deposition at the bottom of a valley (light-shaded square) - 
corresponding to the spin-flip even ^—^i - or particle evap- 
oration from a peak (filled square) - corresponding to i^T- 

In this frame (Fig. ^ lower), a spin flip is equivalent 
to a "particle" which either deposits at the bottom of 
a valley (T^i) or evaporates from a peak (i^T)- In ^ 
single time step each such event occurs with probability 
1/2. For an interface with transverse dimension of order 
L, let us assume that there are of the order of such 
height extrema. Ref. Q predicts /i = 1, but we temporar- 
ily keep the value arbitrary for clarity. Accordingly, in a 
single time step, where all interface update attempts oc- 
cur once on average, the interface center-of-mass moves a 
distance Ay - L'^/^/L. This gives an interface diffusivity 
D ^ (A?/)^ ^ L^^^. We then estimate the lifetime r^iag 
of a diagonal stripe as the time for the interface to move 
a distance of order L to meet another interface. This 
gives Tdiag ~ L'^/D ^ L'^^^. Using the results of Ref. dj, 
we expect Tdiag esc L^. 

The survival probability reflects these two time scales 
(Fig. ^), and their L dependence is cleanly visible in the 
reduced moments = [t^)'^/^ of the time until the 
final state is reached. The main contribution to the mo- 
ments with k < 1 comes from short-lived configurations, 
while for fc > 1 the main contribution comes from long- 
lived diagonal-stripe configurations. Our data for Mk 
with fc < 1 scales approximately as , while for fc > 1, 
Mk scales roughly as L^-^, somewhat faster growth than 
expected (Fig. ||). 

In greater than two dimensions, the probability to 
reach the ground state rapidly vanishes as the system 
size increases. For example, Pgs ~ 0.04 and 0.003 for 



cubic lattices of linear dimension L = IQ and 20. For 
larger lattices, the ground state has not been reached 
in any of our simulations. One obvious reason why the 
system "misses" the ground state is the rapid increase 
in the number of metastable states with spatial dimen- 
sion. This proliferation of metastable states makes it 
more likely that a typical configuration will eventually 
reach one of these states rather than the ground state. 
Another striking feature is that many metastable states 
in three dimensions form connected iso-energy sets, while 
metastable states are all isolated in two dimensions. Thus 
a three-dimensional system can end up wandering forever 
on one of these connected sets. 




FIG. 5. Dependence of = {t'')^^'' on L. Shown are the 
cases fc = 1 (dashed curve), fc = 1/2 (□), 1/4 (A), and 1/10 
(o), as well as fc = 2, 4, and 10 (corresponding filled symbols). 
The thin straight lines have slopes 2 and 3.5. 

A specific example from a simulation on a small cube is 
sketched in Fig. ^. By viewing the spins as cubic blocks, 
the cluster of aligned spins appears as a "building" with 
a 2-storey section (marked 2), an adjacent 3-storey sec- 
tion, and a section (marked oo) which wraps around the 
torus in the vertical direction and rejoins the building on 
the ground floor. The wiggly lines indicate that building 
sections also wrap around in the x- and y-directions. 
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FIG. 6. A stochastic blinker on the cubic lattice. 

The sites marked by x are "blinkers". Consider the 
leftmost such site. Since there are three directions where 
the nearest neighbors are part of the building, this left- 
most spin can flip with zero energy cost. If this occurs, 
then its right neighbor, which was initially stable, can 
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now flip with no energy cost. This sequence can continue 
until the right edge of the oo section of the building but 
no further. Therefore this third-storey addition performs 
a random walk, constrained to move forever in the inter- 
val marked by the x sites. While this construction ap- 
pears idiosyncratic, such a structure is the basic element 
of generic blinkers on even-coordinated Cayley trees. 

Another feature of the final state is that it almost al- 
ways consists of only two interpenetrating clusters which 
both percolate in all three Cartesian directions. These 
two percolating clusters must each contain no convex 
corners to be stable at T = 0. While there are also 
metastable states with many components and metastable 
states with components percolating in one or in two di- 
rections, such configurations are generally not reached 
when the system is large enough. 

The energy decay on the cubic lattice also suggests 
that there exists more than one relaxational time scale. 
Initially, the energy decreases systematically in a man- 
ner consistent with a power-law decay. At longer times, 
however, the energy exhibits plateaux of increasing dura- 
tion, punctuated by small energy decreases. Ultimately 
the final energy is reached, after which constant-energy 
stochastic blinking occurs ad infinitum. Our data for 
the time until the appearance of the first energy plateau 
scales roughly as L^, while the time to reach the final 
energy seems to increase faster than any power of L. 

The phase-space structure of the metastable states ap- 
pears to be a crucial element in understanding the fate 
of Ising ferromagnets. The simplest aspect is to esti- 
mate the number of metastable states Md{N) as a func- 
tion of the spatial dimension d and number of spins 
N ~ U^. In two dimensions, a metastable state con- 
tains alternating horizontal or vertical stripes of up and 
down spins, with each stripe of width > 2. This is 
identical to the number of ground states of a periodic 
Ising chain with nearest-neighbor ferromagnetic inter- 
action Ji and second-neighbor antiferromagnetic inter- 
action J2 (axial next-nearest neighbor Ising (ANNNI) 
model), when J2 = — Ji/2. For the chain with open 
boundaries, the number of metastable state was previ- 
ously found in terms of the Fibonacci numbers For 
the periodic system the number of metastable states is 
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where g = i(l-t-\/5) is the golden ratio, and the first term 
is just the ANNNI model degeneracy. Eq. (|l|) therefore 

gives M2{V) - e^^^ with A2 = ln.g. 

For d — 3, we give a lower bound for the number of 
metastable states by generalizing the stripe states of the 
square lattice. Consider states which consist of an array 



of straight filaments such that each filament cross-section 
is rectangular and the "Manhattan" distance between 
any two rectangles is > 2. The number of packings of 
such filaments is of the order of exp(ci^), where c is a 
constant. This gives the lower bound for the number of 
metastable states, MsiN) > exp(ciV2/3). This same con- 
struction in d dimensions gives Md{N) > exp(ciV^''~^)/''). 
While we have not succeeded in constructing an upper 
bound, it seems plausible that this bound has the same 
functional form as the lower bound; hence Md{N) ^ 
exp {AdN'-'^-^^/'^). In addition, it can be verified that the 
number of metastable states on the Cayley tree grows ex- 
ponentially with the number of spins. Thus metastable 
states become relatively more numerous as the dimension 
increases and their influence on long-time kinetics should 
correspondingly increase. 

In summary, the homogeneous Ising ferromagnet ex- 
hibits surprisingly rich behavior following a quench from 
infinite to zero temperature. On the square lattice, there 
appears to be a non-zero probability of reaching a static 
two-stripe state. Evolution via a diagonal stripe con- 
figuration is responsible for a two-time-scale relaxation 
kinetics. On the cubic lattice, the probability of reaching 
both the ground state or a frozen metastable state van- 
ishes rapidly as the system size increases. The system 
instead reaches a finite iso-energy attractor of metastable 
states upon which it wanders stochastically forever. 
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